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ABSTRACT 


A simple theoretical model is proposed to describe the flow of a turbulent wall jet 
along a curved surface into a quisecent atmosphere. An integral method is used to solve 
the momentum equation and identifies three contributions to the spreading rate of the wall 
jet: (1) turbulent diffusion in the wall jet, (2) wall curvature and (3) rate of change of 

wall curvature. Closed from approximate solutions are found for the case of a plane wall, a 
circular cylinder and a logarithmic spiral surface and comparison with experimental data 
for these three cases is made showing good agreement between theory and experiment. 
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Introduction 


The behavior of a wall jet flowing around a curved surface has been the subject of 
study for almost two hundred years following the first observations of Young* 1 * in 1800. 
Practical application of this phenomena, popularly known as the Coanda effect, have been 
pursued in more recent times with a view to delaying separation of the boundary layer on 
lifting surfaces by using a wall jet over a rounded trailing edge: the wall jet causes the 
flow to remain attached to the surface, displaces the rear stagnation point and induces 
additional circulation, and therefore produces additional lift on the surface. The use of 
this phenomena in the design of circulation control wings (CCW) has received considerable 
attention in recent years and has been described by a number of investigators (see for 
example Wood* 2 * and Englar* 3 *). 

Despite the fact that the phenomena is being used in practical applications to wing and 
rotor design the understanding of the properties of wall jets is still limited and depends 
primarily in experimental information regarding the turbulent mixing, its effect on the 
spreading rate of the wail jet and the corresponding deceleration of the flow as it proceeds 
along the wall. A thorough review was made by Launder and Rodi* 4 * in 1981 of available 
data for wall jets adjacent to surfaces of plane, cylinder and logarithmic spiral surfaces 
and substantial use may be made of this information to distinguish the influence of wall 
curvature, and rate of change of curvature in the stream direction, on the spreading rate 
of the jet. This information together with data from the same review on the mean velocity 
profile in the wall jet provides the basis for checking theoretical models of the flow. 

In the present paper an attempt is made to formulate the simplest possible model for 
the flow of a wall iet emanating from a two dimensional source into quiescent surroundings 
in the presence of a wall of arbitrary shape. The method uses self similar profiles for the 
mean velocity together with a simple eddy viscosity model. The streamwise and radial 
momentum equations are integrated across the wall jet flow to give an expression for the 
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momentum balance including the effect of the pressure gradient induced by the rate of 
change of surface curvature. The streamwise momentum equation is also evaluated at the 
point of maximum velocity to provide a second equation and thereby permit a solution 
for the two unknown quantities b, the jet half width and u TO , the maximum velocity. 
This approach provides approximate closed form solutions for the flow of the wall jet over 
surfaces of various shapes (for quiescent surroundings), and in particular permits a direct 
comparison with the available experimental results for plane, cylindrical or logarithmic 
spiral surfaces. 


Analysis 

The wall jet is considered to comprise two parts: an inner flow adjacent to the wall 
having a highly non-linear velocity profile characteristic of a turbulent wall flow, and an 
outer flow having a velocity profile more typical of a free turbulent plane jet. 

The primary parameters that describe the flow are shown in figure 1. The jet emerges 
from a point source into a fluid at rest and spreads, increasing its width and decreasing 
its velocity due to turbulent diffusion in the jet and friction at the wall. At a distance s 
downstream of the jet exit the velocity to can be expressed as 


u = u m {s)f(y/b) (1) 

where u m is the maximum velocity, occurring at y = y m (s) , and b = 6(s) is the half 
width of the jet (at which point u = lu m ). 

The velocity profile in the outer flow (y > y m ) is assumed to take the form: 


u = u m sech 2 


Ht ~ Sm) 


1~6 


> 6 


( 2 ) 
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where i — yjb. This velocity profile is suggested by the classical free jet solution found 
by Tollmien, modified to give u = u m at £ = The constant k is determined such that 
u = |u m at £ = l(y = b), thus 


k = tanh -1 



= .8814 


(3) 


The velocity profile for the inner flow is assumed to depend on the variable (£/£«,)’* 
as suggested by turbulent wall flow, and is chosen to give a maximum value u = u m at 
the following profile satisfies these conditions 



(4) 


The value £ m , giving the location of maximum velocity, is determined by matching 
the second derivative of the velocity profiles as given by equations (2) and (4), (the first 
derivative is zero since this is the maximum velocity point). 

The result is written 


6.«(1 +*«)-• (5) 

Before proceeding with the analysis, a comparison is made in figure 2 of the velocity 
profiles with the experimental data of Tailland* 5 ** 6 *. The data are typical of those taken 
by a number of authors for wall jet flows over plane and curved surfaces and show that 
the assumed profile given by equations (2) and (4) is reasonable. Additionally in figure 3 
a comparison is made of given by equation (5), with the experimentally determined 
values given by Forthmann* 7 *, Sigalla* 8 * Bradshaw & Gill* s) , Patel (1 °* and Giles* 11 *. Ail the 
data falls within a band .14 < < .16 corresponding to 7 > n > 6 for Reynolds number 
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R, in the range 10 4 to 10 5 ; thus the values of n are consistent with those expected for 
turbulent wall flows. 


In this paper the two primary flow variables to be determined are the jet half width, 
b, and the maximum velocity u m . The continuity and momentum equations for an incom- 
pressible fluid are written 


du d ,, . 

— + = 0 
ds dy ’ 


( 6 ) 


, du , d . —hdp d ,, 2 . 

hu Ts + hV Ty ( ) = TTs + Ty {h r) 


( 7 ) 


and 


u 2 _ 1 ^dp 
R p dy 


( 8 ) 


where 


h = 1 + y/R 


Equations (6) to (8) can be combined to give 


d 2 v?ydR d ,, 2 , 2 d ( py\ 

ds v 1 R 2 ds dy v ds \ p J 


) 


( 9 ) 


neglecting terms of 0 (^r). Integration across the flow gives the integral form of the 
momentum equation: 


d_ 

ds 



1 dR r r°° 2 , 1 / 

+ w77 ll “H = 


4 


( 10 ) 


This equation retains the term which reflects the curvature-induced pressure gradient 
which may be large compared with r„/p, the wall shear stress. 

Substitution of the velocity profiles (equations 2 and 4) into the integrals of equation 
(10) gives 


A n 


-^--j-(6u^) 

u 2 ds v m/ 


+ Cn 


V_dR 
R 2 ds 



(ID 


where 


ji. = ( 1 + |»)(l + *») 1 

C " = \ [(! ‘°8 2 - 1) + 5" + *] + *" 

From equation (11) it can be seen that in the absence of curvature effects the wall jet 
momentum {bu 2 m ) decreases with distance due to skin friction at the wall. However for a 
wall of decreasing radius of curvature (^ < 0) it is possible for the wail jet momentum 
to increase with distance along the wall. This will occur if 


*£(-?) > K 

An approximate form of equation (11) valid for large n and for <7/ — * 0 may be written 


u: 


•* ds K mi k R 2 dS 


(Ha) 


a result which is applicable for large Reynolds number. 

Since there are two unknown quantities in equation (11) it is necessary to use a second 
relationship in order to determine b and u m . It is convenient to use the momentum equation 
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(7) evaluated at the point of maximum velocity, (ie where u = u m , = 0) which may be 
written 


du uv 1 dp 

u i ( — 

^ n 1 1 O 



h dr 

2r 

. m 

.pdy 1 

h pR. 


( 12 ) 

The terms on the left hand side may be written in dimensionless form, neglecting terms 
of 0(b/R) 2 as follows: 


and 


6 ( dn\ 1 

< V ^7/m ” 2 


— —(bu 2 ) - — 
u 2 da ( mJ da 


b /utA _ 1 
n 2 \ R ) u m 


\ 1 

/*»"» 3u , 

) m * 

f/. tM 


1 


2(1 + kn) 


b/R 


±A( 6 U ’)-^ 

u 2 m da V m} da 


b/R 


Finally, 



-t- . A . f-LA (6u * ) _ LIE 1 

1 + kn R [u 2 m da K m) Rda J 


Thus the left hand side of equation (12) becomes 


1 

2 


b/B ) 


1 + kn 


: h-i «) 




-i(* 


1 

1 + kn 



db 

da 
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In the limiting case of large n this expression simplifies to 


1 

2 



db 

ds + 3 k 



In order to evaluate the right hand side of equation (12) it is necessary to express r in 
terms of u through an eddy viscosity e based on a length scale 6(l-f m ) and a characteristic 
velocity u m . It is convenient to write 


€ = ^ p ^ ttm ( 13 ) 

where if is a constant to be determined experimentally. The function gr(^) is chosen such 
that y(£ m ) = 1 and also satisfies the appropriate boundary condition at f = O. 

The shear stress takes the form 



or, substituting for u = u m /(£) and e 


^r = ^(i -Uirg-b/Rfg] (is) 

Similarly, 


b dr 

p<dy 


K 

4*2 




(A + /V) 



(16) 


where / denotes a derivative with respect to 


Since /(£ m ) = ?(&») = 1 and /'(£„,) = O the foregoing expressions evaluated at £ = f m 
combine to give: 


+ =T{^ l ~ (m){l + ub/R)fJl ~ b/Rg ' /f " ) '" + 0WRr } (17) 

The quantity /" is evaluated from 


/ = sech 2 


. i-e . 


as 


aw = 


2k 2 

(1- W a 


Turning nov/ to the function g(£) which describes the variation with £ of the eddy 
viscosity, <?(£) is assumed to take the form: 




( 18 ) 


where the exponent m and the constant a are determined by satisfying the condition 
that 


r ( du u\ 1 _ 

— r- = e [ — ) = — Ct at £ = o 

P u m \dy RJ 2 


Substitution for e and u gives 


K l~£ m 
4* 2 
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at £ = o. 
Thus 



and g’{£m) is evaluated as 


i'(U = (2 - i a + *») 

Thus equation 17 becomes 


( 19 ) 


( 20 ) 


(20a) 


26 


u 


2 

m 



kn 

1 + kn 



l-n/2 + n 2 (l- 
1 - 1 - kn 


b/R 


which has the form, for n large and C/ — » 0; 


26 fhdT_ , 2 T\_ 

< U dy * PR ) ~ 

Substitution into equation (12), and elimination of the quantity jv-^(6u^) using equa- 
tion (11a) gives finally, an expression for the spreading rate — : 



db 

ds 


= K 


+ "M + (i - nr) (*/* §) b > R+0 w R)2 w 

In this equation K is to be determined experimentally for the flat plate (R -* 00 } 
and the term involving (b/R^ j is retained since this may be the same order as ^ (for a 
iog spiral shape, in fact fb /R$) = ^ 
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It is clear from equation 22 that the influence of small curvature (small b/R) may be 
quite significant since the quantity £ b/R appears in the expression for the spreading 
rate. 


Comparison with Experimental Results 

In order to validate the foregoing analysis a comparison is made of the spreading rate 
as given by equation (22) with available experimented results for two cases, namely, the 
circular cylinder (R = constant) and the logarithmic spiral = const ). 

(a) Circular Cylinder 

In this case equation 22 reduces to a linear differential equation for b/R namely 


db f n , , _1 


(22a) 


The solution to this equation may be written in the form 


b 

- = K 


I b/R 


3 log(l + f b/R) 


or, neglecting terms of 0(b/R) 2 , 


(23) 


i/3 = *( 1 + s 6/R ) 


(23a) 


The constant K is determined from the flat plate wall jet (b/R = O) as K = .07, and 
equation (23a) with k = .8814 gives 


b/s = 07(1 + 4b/ R) (n = 7) 


(236) 
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6/s = .07(1 + 3.46/5) (n = 6) 


(23c) 


These expressions are shown in figure 4 together with the available experimental results 
of Fekete* 12 * , and Wilson and Goldstein* 13 ). In this regard the review by Launder and 
Rodi* 14 ) suggest that the data of Fekete is possibly more representative of two dimensional 
flow than that of Wilson and Goldstein. It is in fairly good agreement with the theoretical 
expression for n = 6 (the Reynolds number is in the range 4 to 13 x 10 s ). 

(b) Logarithmic Spiral 

For a logarithmic spiral R/s = constant and equation (22) has a self similar solution 


6/s = 


k [l + ^ 6/5 ] 
1 - ii^6/5 


» .07 


[1 + 86/5j 
[1 - 6/5] 


(n = 7) 


.07(1 + 6.86/5 ] 

[I-*/*] 


(n = 6) 


(24) 


(24a) 


(246) 


A comparison of these curves with the experimental results of Guitton and Newman* 14 ) 
and Kamemoto* 15 ) is made in fig. 5. Here Guitton and Newman’s results are considered to 
be most representative of two dimensional flow. Equation (24) reflects the much greater 
increase of b/s with ’b/R for the logarithmic spiral, compared to the circular cylinder 
equation (23a), and this trend is also seen in the experimental data. 


Conclusions 

From the foregoing analysis and comparison with experimental results it is seen that 
the wall jet is influenced by three effects: 
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(a) The turbulent diffusion in the wall jet giving rise to a linear spreading rate dbjds = 
K « .07. 

(b) A coupling between the eddy viscosity and the curvature of the wall arising from the 
term eu/R ) in the stress gradient and yielding a contribution of = K^b/R. 
(typically this contribution is 0(^b/R) for Reynolds numbers typical of laboratory 
experiments), and 

(c) A curvature-induced streamwise pressure gradient giving rise to a contribution of 
0(b/R ) 2 ^ 7 - When ji < O this effect causes the jet momentum ( 6u £,) to increase 
in the streamwise direction and help confine the jet to a thin layer as it proceeds along 
the surface. This is an inertial term, largely independent of Reynolds number, and is 
an essential feature of Coanda turning of the flow. 

These influences are additive and give an approximate relationship between the spread- 
ing rate £, the half width of the jet and the surface radius of curvature R, of the form 


db 

ds 


K 0 + K l - + K 2 

it 



dR 

ds 


which can be integrated to give b(s) when the wall shape R(s) is known. 

The use of this result for the circular cylinder (R = const) and for logarithmic spiral 
surfaces (b/R = const) shows good agreement with the previous experimental results of 
several authors. This provides some confidence that a simple algebraic eddy viscosity 
model, used in conjunction with the appropriate expression for the shear stress (including 
curvature terms), is sufficient to describe the primary features of the wall jet. 
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as a function of b/R for circular cylinders. 
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Figure 5. Comparison of theory and experiment for the wall jet spreading rate b/S 
as a function of b/R for logarithmic sprial surfaces. 


